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We study the complex dynamics of a one-dimensional Bose gas subjected to a dissipative local 
defect which induces one-body atom losses. In experiments these atom losses occur for example when 
a focused electron or light beam or a single trapped ion is brought into contact with a quantum gas. 
We discuss how within such setups one can measure or manipulate densities locally and specify the 
excitations that are induced by the defect. In certain situations the defect can be used to generate 
entanglement in a controlled way despite its dissipative nature. The careful examination of the 
interplay between hole excitations and the collapse of the wave function due to non-detection of loss 
is crucial for the understanding of the dynamics we observe. 



Ultracold quantum gases offer the possibility to ad- 
dress many open questions from various areas of physics 
due to their excellent control and fast tunability. 
Achievements of the last decade range from the observa- 
tion of strongly correlated one-dimensional bosonic gases, 
so-called Tonks-Girardeau gases [l|, to the realization 
of a Mott-insulating state in optical lattices, both for 
bosonic and fermionic Q atoms. Only very recently, 
new experimental setups opened the way to resolve sin- 
gle atoms in quantum gases. Fluorescence techniques 
4 , [ij or scanning with a highly focused electron beam 
6| can address single columns in the quantum gas with 
a resolution of the optical lattice spacing. Real three- 
dimensional resolution could be reached using a trapped 
ion 0, Q . These local techniques can be employed to de- 
tect or even to manipulate quantum gases. Consequently, 
these tools enable one to take advantage of the presence 
of the trap to access a large range of the homogeneous 
phase diagram in a single realization of an experiment. 

The main subject of this work is to explore how a 
highly localized dissipative defect that generates atom 
losses, such as an ion or a focused electron or light beam, 
can be used in a one-dimensional system to probe and 
to manipulate properties of the quantum state. For in- 
stance, we show how to measure the local density effi- 
ciently and explain how shock waves emerge in superfluid 
states. At strong dissipative strength subsystems left and 
right from the defect separate spatially but entanglement 
between the subsystems can be manipulated in a con- 
trolled way and is directly related to the experimentally 
measurable atom loss. Our theoretical model is a cor- 
related onc-dimcnsional bosonic lattice system which is 
coupled to a Markovian environment by one-body atom 
losses. Actually, the situation of a local defect immersed 
in a quantum system is very general and not only of in- 
terest in the context of cold atomic gases but also in con- 
densed matter physics. It represents a non-trivial many- 
body problem and in order to capture the coaction of 
continuously created losses and intrinsic dynamics of the 
correlated quantum states correctly, the full treatment 



of the master equation of the interacting system beyond 
mean- field 13, Hi is crucial [l2j . 

In a broad parameter regime, a one-dimensional Bose- 
gas which is subjected to an optical lattice potential can 
be well described by the Bose-Hubbard model, 

H = J2-J Ube+i + bl+^be) + ^ni{n, - 1). (1) 



Here b\ {bg) are the bosonic creation (annihilation) oper- 
ators at site £ and is the density operator. We assume 
a system with open boundaries; integer indices £ ranging 
from — (L — l)/2to (L — l)/2 represent individual sites 
that are spaced by distance a. The first term models 
the kinetic energy of the atoms in the periodic poten- 
tial and the second term the on-site interaction due to 
s-wave scattering. Initially, in the absence of a dissipa- 
tive defect, the atomic cloud is prepared in its ground 
state. For weak interaction this ground state is a super- 
fluid with large on-site density fluctuations. At integer 
filling the state is Mott-insulating for interactions above 
a critical strength (in ID {U/J)c ~ 3.4 at filling n = 1). 
In contrast, at non-integer filling the weakly interacting 
superfluid is connected by a crossover to the strongly in- 
teracting Tonks-Girardeau gas of impenetrable bosons. 

At time t = Q the dissipative defect (localized at site 
t — Q), e.g. a trapped ion or an electron or light beam, 
is brought into contact with the bosonic cloud and starts 
generating atom losses at the central lattice site only. 
Experiments have reached or are close to reach such res- 
olution 0, S, 0] ■ From the theoretical point of view it is 
of secondary importance if one or few sites are affected. 
We assume that atoms are expelled quickly (compared 
to the timescales of the quantum gas) to the free space 
continuum and that the defect itself does not become 
correlated to the quantum gas. For this Marko- 
vian approximation can be used to derive a master equa- 
tion [13| for the density matrix p of the bosonic cloud, 
p = Cf' p -I- CP p., Cf' p := —i\H,p\, with the dissipator 
CP p := Tbopb^ — ^710/9 — ^priQ. The strength of the cou- 
pling to the environment F depends on the cross-section 
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of the scattering processes between the atoms and the de- 
fect. In the ion-atom hybrid system the loss mechanism 
can be seen as an elastic scattering event of an atom with 
a hot ion. In principle, its strength can be widely tuned 
by using a suitable magnetic Feshbach resonance [13] or 
using different kinds of atoms and ions. Typical ener- 
gies of an ion are of the order of Ifc^K, while the atomic 
trapping potential is of the order of Ifc^/xK There- 
fore the atoms will be scattered to the free continuum 
in a time much shorter than the many-body timescales 
in the cloud (e.g. J ^ 10~*fcBK). When shining a fo- 
cused electron beam on the cold atomic sample [6[ the 
loss mechanism corresponds to an ionization process, for 
which a similar separation of energy scales is present (in 
this case F can be tuned by the intensity of the beam). 
Also in experiments with fluorescent light Q the energy 
gain by spontaneous emission is sufficient to transfer an 
atom to the continuum or at least to higher bands. Con- 
trarily to previous experiments Q in which the lattice 
was increased prior to illumination in order to avoid this 
process, in the current proposal the lattice is kept con- 
stant and the focus is on the induced dynamics. 

The non-equilibrium problem of the master equation 
is solved by sampling over quantum trajectories 13, is} . 
The propagation of the stochastic wave functions with 
the non-hermitian Hamiltonian H^s = H ~ iVrio [l5j is 
achieved by the time-dependent density matrix renormal- 
ization group (DMRG) The imaginary part of 

Hes represents the effect of non-detection of atom loss. 
Detection of atom loss is simulated by annihilation of 
bosonic fields at stochastically chosen times. For the 
study of the dynamics of locally created defects in one- 
dimensional systems, DMRG works quite efficiently as 
seen in a spin system in Ref. [I8j . In the present context, 
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bond dimensions 100-200 in the DMRG and time-steps 
of 0.001 — 0.005 h/J lead to errors negligible compared 
to the statistical error of the sampling procedure. The 
optimal choice of DMRG parameters is important, since 
averaging over more than 6000 trajectories may be neces- 
sary. Quantum trajectories allow for a natural definition 
of entanglement entropy in this open system fioj . One 
defines Si , the entropy of entanglement for a bipartition 
of the system at a site I, by the average over the von 
Neumann entanglement entropies of single trajectories. 

One of the experimentally accessible quantities is the 
total atomic loss N{t) ~ \ni^') — ng(t = 0)] measured 
after a given time t (n((t) ~ {ng){t)). On short time 
scales, as long as current flows into the central site are 
negligible (guaranteed for t h/ J, t ^/r), the total 
loss is 



N{t) 



.(1 



where m, 



no{t = 0). (2) 



We will focus on the response of the system for times be- 
yond this short time behaviour in which a complex many- 
body dynamics sets in that depends on the nature of the 
underlying quantum state in a non-obvious manner. 
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FIG. 1. (Color online) Weak dissipative coupling T/J — 0.25, 
L = 105 and 53 or 105 particles. In (a) and (b) we show 
density profiles at different times. Densities are rescaled by 
their initial values and size-dependent features in the vicinity 
of the edges of the systems are discarded. We depict the 'light- 
cone' of waves moving with the sound velocity [2y] by vertical 
lines. Star and cross symbols represent |i/)nd(i)), IV'hoiG(i)) 
respectively. In the lower panels we plot the total atom loss 
(c) and the central density (d), both rescaled by the initial 
density nin. In (c) 4 curves almost lie on top of each other. 
Boundary effects are eliminated by interrupting simulations 
before reaching the recurrence times. Statistical errors are 
either marked by bars or smaller than line width or symbol 
size. 



First we analyse the time evolution of the density pro- 
files in the case of a weak dissipative defect, T / J = 0.25 
[Fig. [TJa-b)]. For the strongly interacting Tonks gas 
[J7/J = 8, = 0.5 in Fig.[T](b)] a dark density modula- 
tion moving at the speed of sound (close to 2Ja/h [2^) is 
the dominating feature in the relaxation. Nearly identical 
evolutions of the density were found for a Mott insulat- 
ing state (not shown), in which hole modes [2l| are also 
propagating with velocity v « 2Ja/h. The insulating na- 
ture does not play an important role for the transport 
of the excitation. In contrast to these strongly inter- 
acting cases, in the weakly interacting superfiuid state 
\U / J = 1, = 0.5 in Fig. [Ija)] we observe the forma- 
tion of a remarkably strong bright shock wave at a speed 
which strongly exceeds the sound velocity. Only after this 
density wave the expected dark modulation propagates. 

To reveal the origin of this shock wave we analyze in 
Fig- [Ha) density profiles of single trajectories which allow 
to separate the effects of the detection and non-detection 
of atom losses: (i) the wave function evolves in time by 
the Hamiltonian and the non-detection part of the dis- 
sipative action only, |V'nd(i)) = e^'^^'="*^^\'ipo); (ii) the 
wave function is the evolution of a hole, generated at 
t = by a single local atom loss, with the Bose-Hubbard 
Hamiltoninan H]), |i/'hoic(t)) = e"*^*/'^6o|V'o)- From this 
analysis follows clearly that the continuously induced col- 
lapse of the wave function by the non-detection of atom 
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losses is the driving force behind the shock waves. The 
propagation of a locaUzed hole does not lead to signifi- 
cant shock wave formation although in the present case 
the dark perturbation in the density of |V'hoic(*)) is much 
larger than in \tpnd{t))- The highly non-linear effect of 
the non-detection of loss is directly related to the non- 
local correlations present in the superfluid and therefore 
not observed in the more strongly interacting cases. 

Surprisingly almost no sign of this very different dy- 
namics due to the underlying quantum phase can be no- 
ticed in the evolution of the total atom losses [Fig. [IJc)]. 
At short times, before a current sets in, N{t) should obey 
the exponential law ([2]) for any type of system, which is 
indeed the case up to t « h/AJ. More interestingly, at 
larger times N{t)/nin shows a linear increase with prac- 
tically no dependence on the initial state. In Fig. (Hd) 
the central density no{t) is plotted in order to visualize 
the tiny differences in the behaviour of the total loss (it 
represents the loss rate by the relation N{t) = rno(i)). 
no{t) describes how, through a slowly decaying, oscilla- 
tory transient regime, a quasi-steady value is reached. 
In order to investigate systematically the onset of cur- 
rent we extract a stationary value no by averaging over 
the time interval lOh/J < t < 15h/J. Looking at no 
for various dissipative coupling strength [Fig. EJa)] we 
found a generically obeyed relation no/nin ~ e~'"/^'^ up 
to coupling strength F/J w 0.4. From this result we can 
conclude that H/2J is the characteristic timescale for the 
onset of a current. From our analysis of the density pro- 
files we can understand this time scale from the modes 
that dominate the relaxation process. Indeed, from the 
velocities of sound and 'hole' modes observed in Tonks 
and Mott systems {v w 2Ja/h) a relaxation time h/2J 
would be expected. In principle, for the weakly interact- 
ing superfluid one has a much smaller sound velocity 
{v ~ l.lJa/h in the case plotted in Fig. [BJa)). However, 
we found that the faster shock-waves are the dominant 
modes, which apparently again lead to a relaxation time 
scale h/2J. 
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FIG. 2. (Color online) (a) no in the weakly dissipative regime. 
Error bars indicate systematic fluctuations and statistical er- 
rors, (b) Long-time loss rate for the strongly dissipative case, 
(c) Loss rate from weak to strong dissipation. Lines are guides 
to the eye. 



Hence, the short-time exponential decay ([2|) followed 
by a long time linear rise with slope ^ e^^^'^ is generic 
and obeyed for very different types of systems. This 
means that one can actually infer the value of the initial 
density from atom loss measurements without explicit 
knowledge of the nature of the state of the sample be- 
yond the short time regime ([2]). Long exposure of the 
sample to the defect is helpful to increase the detected 
signal. 
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FIG. 3. (Color online) Numerical results for strongly dissi- 
pative coupling, r/J — 8. For explanation of (a-c) see Fig. 
[T] Panel (d) shows the evolution of the von-Neuman entropy 
between the subsystems separated by the defect. 

Let us turn to the response of the atomic cloud to a 
strong dissipative defect (F :S> J,U). The density pro- 
files [Fig. [2Ia-b)] show the propagation of modes that 
are very similar to the weakly dissipative case, but with 
more pronounced density modulations. At the central 
site the strong dissipation leads to a very rapid suppres- 
sion of the density suggesting a decoupling of the system 
at the defect. In contrast to the previously considered 
weakly dissipative case, the atom loss is now highly sen- 
sitive on the initial state [Fig. [3][c)]. The von Neumann 
entropy Sq also strongly depends on the initial state and 
shows quite unexpected behaviour: strong entanglement 
between the two separated subsystems can be generated 
by the defect [Fig.|3Id)]. Only at short times the entropy 
of entanglement is reduced by the projection performed 
by the measurement. At larger times the entanglement 
grows approximately linearly in time. In the strongly in- 
teracting systems the linear rise is found numerically to 
be proportional to the atom loss S'o(t) ~ aN{t), a < 0.5, 
where a maximal efficiency a = 0.5 is reached in the 
hard-core limit J7/J — > oo (not shown). For the curves 
at strong but finite interactions plotted in Fig. [Sjd) we 
find a « 0.3. The effect of non-detection of atom loss 
is crucial for this specific relation between loss and en- 
tanglement entropy. If there was only creation of holes 
at different instances of times we would expect that each 
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hole, which splits into two mutually entangled excitations 
propagating into opposite directions, would lead to an in- 
crease in the entanglement entropy by one unit. However, 
the collapse of the wave function by non-detection of loss 
destroys such pairs and apparently limits the efficiency 
of entanglement generation to a < 0.5. When decreasing 
the interaction, efficiency a is lowered and the production 
of entanglement becomes nearly completely suppressed 
in the weakly interacting superfluids [Fig. Eljd)]. This 
could be explained by the presence of density fluctuations 
in the initial wave function. The propagating hole pairs 
do not create new highly excited states, as it is the case 
at strong interactions, but only leads to redistribution of 
weights between states. 

It is interesting to gain an analytical understanding of 
the long time-limit of the density profiles and the corre- 
sponding atom loss, especially because of their relation to 
the entanglement entropy in strongly interacting states. 
A perturbative analysis by adiabatic elimination of oc- 
cupations of the central site (see e.g. [lij) can be per- 
formed. From the effective master equation so derived, 
N{t -^h/J) = 2^,^±l2n,(t) + (6t6_,)(t)]+ O (^), it 

follows that the loss rate decreases at increasing coupling 
^9 — — 

This is a counterintuitive result 
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strength as 

which is known for discrete measurements as the quan- 
tum Zeno effect [23|. Using that the density close to the 



center can be described by a static boundary situation 
and correlations between the two sites next to the de- 
fect can be neglected, the perturbative expression can be 
further simplified, 

8 P 

N{t:^h/J)^—n±^L-i)/2{t = 0). (3) 

A comparison between this expression and the numerical 
results gives an excellent agreement down to a dissipa- 
tive coupling of F w 8J for strong interactions [Fig.[2][b)] 
(the loss rate is extracted by averaging over the station- 
ary state for lOh/J < t < Ibh/J.). In contrast, the 
agreement to the weakly interacting superfluid is not 
yet fully established. The averaged values nevertheless 
obey the characteristic J^/F scaling with a strong ten- 
dency towards the stationary rate. We expect that for 
longer times than were accessible to us, this value will be 
reached. 

In conclusion, we have shown that dissipative coupling 
can be used to measure densities accurately and to ma- 
nipulate properties of a quantum gas in a controlled way, 
although it induces highly non-trivial many-body dynam- 
ics. Rather unexpected is the strong growth of entangle- 
ment between seemingly decoupled subsystems. It is very 
interesting for practical applications that this entangle- 
ment can be manipulated, and, at the same time, be mon- 
itored by the experimentally easily accessible atom loss. 
We have calculated these effects in a one-dimensional sys- 
tem with great precision, but the main features, e.g. den- 
sity measurements or the Zeno-efFect, are expected to 



occur in higher dimensions as well. Further, entangle- 
ment generation should be possible analogously when 
placing a one-dimensional defect into a two-dimensional 
lattice. In experiments, one interesting avenue to follow 
now would be to use a moving defect to measure static 
or time-resolved correlation functions or to manipulate 
different regions of the system. Moreover, the loss mech- 
anism can be used to remove entropy-rich regions and 
thus may serve as a tool to further cool down an atomic 
gas. The new experimental tools that manipulate quan- 
tum gases by local dissipation open a new perspective on 
many-body dynamics of correlated systems. 
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